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Abstract 


Recently, a new framework for describing the multiverse has been proposed which is based 
on the principles of quantum mechanics. The framework allows for well-defined predictions, 
both regarding global properties of the universe and outcomes of particular experiments, ac- 
cording to a single probability formula. This provides complete unification of the eternally 
inflating multiverse and many worlds in quantum mechanics. In this paper we elucidate how 
cosmological parameters can be calculated in this framework, and study the probability dis- 
tribution for the value of the cosmological constant. We consider both positive and negative 
values, and find that the observed value is consistent with the calculated distribution at an 
order of magnitude level. In particular, in contrast to the case of earlier measure proposals, our 
framework prefers a positive cosmological constant over a negative one. These results depend 
only moderately on how we model galaxy formation and life evolution therein. 


1 Introduction 


An explanation of a small but nonzero cosmological constant is one of the major successes of the 
picture that our universe is one of the many different universes in which low energy physical laws 
take different forms [I]. Such a picture is also suggested theoretically by eternal inflation [2] and the 
string landscape [3]. This elegant picture, however, has been suffering from the predictivity crisis 
caused by an infinite number of events occurring in eternally inflating spacetime. To make physical 
predictions, we need to deal with these infinities and define physically sensible probabilities [4]. 

Recently, a well-defined framework to describe the eternally inflating multiverse has been pro- 
posed based on the principles of quantum mechanics [5]. In this framework, the multiverse is 
described as quantum branching processes viewed from a single “observer” (geodesic), and the 
probabilities are given by a simple Born-like rule applied to the quantum state describing the entire 
multiverse. Any physical questions—either regarding global properties of the universe or outcomes 
of particular experiments—can be answered by using this single probability formula, providing 
complete unification of the eternally inflating multiverse and many worlds in quantum mechanics. 
Moreover, the state describing the multiverse is defined on the “observer’s” past light cones bounded 
by (stretched) apparent horizons; namely, consistent description of the entire multiverse is obtained 
in these limited spacetime regions. This leads to a dramatic change of views on spacetime and 
gravity. 

In this paper we present a calculation of the probability distribution of the cosmological con- 
stant in this new framework of the quantum multiversel{ We fix other physical parameters and 
ask what values of the cosmological constant A we are likely to observe. In Section [2] we begin 
by reviewing the proposal of Ref. [5], and we then explain how cosmological parameters can be 
calculated in Section 3] While the framework itself is well-defined, any practical calculation is nec- 
essarily approximate, since we need to model “experimenters” who actually make observations. In 
our context, we need to consider galaxy formation and life evolution therein, which we will do in 
Section [4] We present the result of our calculation in Section [5] We find that, in contrast to the 
case with some earlier measures [li], the measure of Ref. does not lead to unwanted preference 
for a negative cosmological constant—in fact, a positive value is preferred. We find that a simple 
anthropic condition based on metallicity of stars is sufficient to make the calculated distribution 
consistent with the observed value at an order of magnitude level. We conclude in Section [6] 

Appendix [A] lists formulae for galaxy formation used in our analysis. Appendix [B] discusses the 
anthropic condition coming from metallicity of stars. 


'For earlier studies of the cosmological constant in the context of geometric cutoff measures, see Refs. [6]—[IQ]. 


2 Probabilities in the Quantum Mechanical Multiverse 


Here we review aspects of the framework of Ref. [5] which are relevant to our calculation. In 
this framework, the entire multiverse is described as a single quantum state as viewed from a 
single “observer” (geodesic). It allows us to make well-defined predictions in the multiverse (both 
cosmological and terrestrial), based on the principles of quantum mechanics. 

Let us begin by considering a scattering process in usual (non-gravitational) quantum field 
theory. Suppose we collide an electron and a positron, with well-defined momenta and spins: |ete~) 
at t = —oo. According to the laws of quantum mechanics, the evolution of the state is deterministic. 
In a relativistic regime, however, this evolution does not preserve the particle number or species, so 
we find 

W(t =—-00)=|ete™) + Wt=+00) =e leTe )+ ey lutu)t+---, (1) 


when we expand the state in terms of the free theory states (which is appropriate for t + --oo when 
interactions are weak). The Hilbert space of the theory is (isomorphic to) the Fock space 


H= Oxi, (2) 
n=0 
where Hp is the single-particle Hilbert space. Various “final states,” |ete~) , |u* 7) ,---, in Eq. 


arise simply because the time evolution operator causes “hopping” between different components 
of the Fock space in Eq. (2). 

The situation in the multiverse is quite analogous. Suppose the universe was in an eternally 
inflating (quasi-de Sitter) phase © at some early time t = to. In general, the evolution of this state 
is not along the axes determined by operators local in spacetime. Therefore, at late times, the state 
is a superposition of different “states” 


Wit=to)=|4) 7 Wth= Silt) |(cosmic) configuration 7) , (3) 
when expanded in terms of the states corresponding to definite semi-classical configurations. The 
Hilbert space of the theory is (isomorphic to) 


H= ‘) TLAas Hy = FLA walle & Th heen (4) 
M 


where Hy, is the Hilbert space for a fixed semi-classical spacetime M, and consists of the bulk 
and horizon parts Hwbux and Hm horizon. (The quantum states are defined on the “observer’s” 
past light cones bounded by apparent horizons.) The final state of Eq. becomes a superposition 
of different semi-classical configurations because the evolution operator for V(t) allows “hopping” 
between different Hy, in Eq. (4). 


As discussed in detail in Ref. [5], any physical question can be phrased as: “Given what we know 
about our past light cone, A, what is the probability of that light cone to have properties B as 
well?” This probability is given by 

dt (U(t)| Oane |W (t)) 
ra¢ PV aaea ee tae 6) 
Jat (UO)| Oa |¥@) 
assuming that the multiverse is in a pure state |W(t)). (The mixed state case can be treated 
similarly.) Here, O,4 is the projection operator 


Os = Slog) (aul (6) 
where |a4;) represents a set of orthonormal states in the Hilbert space of Eq. (4), i.e. possible past 
light cones, that satisfy condition A (and similarly for O,qg). Despite the fact that the t integrals 
in Eq. () run from t = to to ov, the resulting P(B|A) is well-defined, since |W(t)) is continually 
“diluted” into supersymmetric Minkowski states [5]. 

The formula in Eq. (or its mixed state version) can be used to answer questions both regarding 
global properties of the universe and outcomes of particular experiments. This, therefore, provides 
complete unification of_the two concepts: the eternally inflating multiverse and many worlds in 
quantum mechanics To predict /postdict physical parameters x, we need to choose A to select 
the situation for “premeasurement” without conditioning on z. We can then use various different 
values (ranges) of x for B, to obtain the probability distribution P(x). In the next section, we 
discuss this procedure in more detail, in the context of calculating the probability distribution of 
the vacuum energy, « = px = A/87Gy. 


3 Predicting/Postdicting Cosmological Parameters 


In order to use Eq. to predict/postdict physical parameters, we need to know the relevant 
properties of both the state |W(t)) (or its bulk part Ppux = TYhorizon |W(t)) (Y(t)]) and the operators 
Ox, and Oxynpg. Here we discuss them in turn. 

In general, the state |W(t)) depends on the dynamics of the multiverse, including the scalar 
potential in the landscape, as well as the initial condition, e.g. at t = tg. Given limited current 
theoretical technology, this introduces uncertainties in predicting physical parameters. However, 
there are certain cases in which these uncertainties are under control. Consider x = py. We are 

?The claim that the multiverse and many worlds are the same has also been made recently in Ref. [12], but the 
physical picture there is very different. Those authors argue that quantum mechanics has operational meaning only 
under the existence of causal horizons because making probabilistic predictions requires decoherence with degrees of 


freedom outside the horizons. Our picture does not require such an extra agent to define probabilities (or quantum 
mechanics). The evolution of our W(t) is deterministic and unitary. 


interested only in a range a few orders of magnitude around pa op, & (0.0024 eV)* [3], which is 
tiny compared with the theoretically expected range —Mf, S pa S M5). Therefore, unless the 
multiverse dynamics or initial condition has a special correlation with the value of the vacuum 
energy in the standard model (SM) vacua, we expect that the probabilities of having these vacua in 
|W(t)) is statistically uniform in 2 within the range of interest. (This corresponds to the standard 
assumption of statistical uniformity of the prior distribution of p, [i].) The distribution of x = pa is 
then determined purely by the dynamics inside the SM universes, i.e., the probability of developing 
experimenters who actually make observations of the vacuum energies. 

Let us now turn to the operators O, and Oynqg. In order to predict the value of the vac- 
uum energy which a given experimenter will observe, we need to choose Oy to select a particular 


“premeasurement” situation for that experimenter, i.e. 


A: aparticular “premeasurement” situation 


P(pn) don = P(BIA), { B: py, < vacuum energy < pa + dpa, (7) 


where P(B|A) is defined in Eq. (6). Here, we have assumed that the number of SM vacua is suffi- 
ciently large for p, to be treated as continuous in the range of interest. In general, the specification 
of the premeasurement situation can be arbitrarily precise; for example, we can consider a particular 
person taking a particular posture in a particular room, with the tip of the light cone used to define 
|W(t)) located at a particular point in space. In practice, however, we are interested in the vacuum 
energy “a generic observer” will measure. We therefore need to relax the condition we impose as 
A; in other words, we need to “coarse grain” the premeasurement situation. In fact, some coarse 
graining is always necessary when we apply the formalism to postdiction (see discussions in Ref. [5]). 

What condition A should we impose then? To address this issue, let us take the semi-classical 
picture of the framework, discussed in Section 2 of Ref. [5]. In this picture, the probability is given 
by 


ag ae 


where M4 is the number of past light cones that satisfy A and are encountered by one of the n 
geodesics emanating from randomly distributed points on the initial hypersurface at t = to. (This is 
equivalent to Eq. in the regime where the semi-classical picture is valid.) Since we vary only pa, 
all the SM universes look essentially identical at early times when the vacuum energy is negligible. 
The assumed lack of statistical correlation between py, and the multiverse dynamics then implies 
that we can consider a fixed number of geodesics emanating from a fixed physical volume at an 
early time (e.g. at reheating) in these universes, and see what fraction of these geodesics find the 
“premeasurement” situation A in each of these universes. 

Given that we are focusing on the SM universes in which only the values of p,q are different, it 
is reasonable to expect that all the experimenters look essentially identical for different p,, at least 


statistically—in particular, we assume that they have similar sizes, masses, and lifetimes. With this 
“coarse graining,” the condition A can be taken, e.g., as: the geodesic intersects with the body of 
an experimenter at some time during their life. In practice, this makes the probability proportional 
to the fraction of a fixed comoving volume at an early time that later intersects with the body of 
an observer. Note that the details of the condition A here do not matter for the final results—for 
example, we can replace the “body” by “head” or “nose” without changing the results because its 
effect drops out from the normalized probability. Thus, in this situation (and any situation in which 
condition A can be formulated entirely in terms of things directly encountered by the geodesic), 
the semi-classical approximation to the scheme of Ref. can be calculated as the fat geodesic 
measure outlined in Ref. 3] We emphasize that the consistent quantum mechanical solution to 
the measure problem in Ref. [5] forces this choice on us. 

We can now present the formula for P(p,) in a more manageable form. Since the probability 
for one of the geodesics to intersect an experimenter is proportional to the number of experimenters 
and the density of geodesics, we have 


P(py) x > obec Peeod,a; (9) 


a€ habitable galaxies 


where Nops,a and Pgeod,a are, respectively, the total number of observers/experimenters and the 
density of geodesics in a “habitable” galaxy a. Here, we have approximated that Pgeoa,a is constant 
throughout the galaxy a. Note that since we count intersections of experimenters with geodesics, 
rather than just the number of observers (as in much previous work, e.g. [6]), our results differ from 
such previous results by our factor of Pgeoaa. Our remaining task, then, is to come up with a scheme 
that can “model” Nobs,¢ and Pgeod,a reasonably well so that the final result is not far from the truth. 


4 Approximating Observers 


In this section, we convert Eq. (9) into an analytic expression that allows us to compute P(p,) 
numerically. We focus on presenting the basic logic behind our arguments. Detailed forms of the 
functions appearing below, e.g. F(M,t) and H(t’; M,t), as well as useful fitting functions, are given 
in Appendix [A] 

Let us begin with Nopsa. We assume that, at a given time t, the number of observers arising in 
a given galaxy a is proportional to the total number of baryons in a 


CN eked 
dt 


3To our knowledge, no detailed study of the probability distribution of the cosmological constant according to the 
fat geodesic measure has been published prior to this work. 


(t) & Nz,(t), (10) 


as long as stars are luminous. This assumption is reasonable if the rate of forming observers is 
sufficiently small, which appears to be the case in our universe. To estimate the number of baryons 
existing in galaxies, we use the standard Press-Schechter formalism [15], which provides the fraction 
of matter collapsed into halos of mass larger than M by time t, F(V/,t). Since the amount of baryons 


collapsed is proportional to that of matter, we can use this function F’ to estimate the number of 


zal M, t) 
as far i im ) peeoa(M, t). (11) 


The expression of Eq. (LL) does not take into account the fact that forming intelligent observers 


observers and find* 


takes time, or that a appear only when stars are luminous (which we postulate, motivated 
by the assumption that we are typical observers). To include these effects, we use the extended 
Press-Schechter formalism [16], which gives the probability H(t’; M,t) that a halo of mass M at 
time t virialized before t’. The probability density P(p,) can then be written as 


dF(M,t) 
s— fu fav ) fit ose WE t) — H(t — thum; M, t)} Peeoa(M, t), (12) 


where teyo, and tpurm are the time needed for intelligent observers to evolve and the characteristic 
lifetime of stars which limits the existence of life, respectively. 

The density of geodesics pgeoa(M,t) is proportional to that of a dark matter halo of mass M at 
time t, which is given by its average virial density: 


dF(M,t.)\ [*., » @2F(M,t’) 
x | ————_ : ——— 1 
Pecoa(M, t) ( dM ) / dt pat) Tra oy 


where ¢t, = min{t, tstop(Z)} with tstop(JZ) the time after which the number of halos of mass M 
starts decreasing, i.e. when merging into larger structures dominates over formation of new halos: 
d?F/dMdt\t=t,.)(@) = 0. (For the explicit expression of pyir, see Appendix [A]) In the interest of 
speeding up numerical calculation, we approximate this by the virial density at the time when the 
rate of matter collapsing into a halo of mass M, i.e. —d?F'/dMdt, becomes maximum: 


Peeoa(M, t) x Pvir(T(M)) ; (14) 
where 7(M) is given by 
3F(M,t 
a F(M, t) Zi, (15) 
dMdt? |,_,(a 


This approximation is indeed reasonable at the level of precision we are interested in: it works at 
the level of 20% for t 2 1.77(M) where the contribution to P(p,) almost entirely comes from. 


4Note that the sign of dF/dM is negative because of the definition of F. 


Finally, there will be several additional anthropic conditions for a halo to be able to host intelli- 
gent observers. For example, the mass of a halo may have to be larger than some critical value Min 
to efficiently form stars [17], and smaller than M,,,x for the galaxy to be cooled efficiently [79]. 
Considering these factors, we finally obtain from Eqs. and (14) 


1 ft. fe dr Mt 
P(p,) = -— | at f ay AOE) {H(t — tovor; M,t) — H(t — thurs M, t)} pyir(t(M)) n(M, 2), 
t 


evol Arnin 
(16) 
where JN is the normalization factor. Here, 


oe) for px => 0 
th = 


_ f[ 17 
teeanch a 6GylPal for PA < 0, ( ) 


and we have put anthropic conditions besides Miinmax in the form of a function n. Note that F, 
H, and pyir (and possibly n) all depend on the value of the vacuum energy pa; see Appendix [A] 

In summary, (dF'/dM) (A |+-+4,,.. — Hlt-n,,,,,) 7 dM dt counts the (expected) number of observers 
in halos with mass between VM and M+dM at time between ¢ and t+dt, and pyir(7) is proportional 
to the density of geodesics in such a halo and time, and so Eq. (16) gives the probability by 
counting their intersections (as in Eq. (Q)), with n implementing some anthropic conditions. One 
well-motivated origin for n is metallicity of stars, which affects the rate of planet formation (see 
e.g. Refs. [21]). Here we simply model this effect by multiplying some power m of integrated 
star formation up to time t — teyo1, which we assume to be proportional to the integrated galaxy 
formation rate for M > Minin: 


n(M,t) & (F(Mnin, min{t — tevol, tps) — F(M, min{t — tev, baton} )) 3 (18) 


where tstop is determined by d{ F(Mmin, t') — F(M, t’)}/dt'|y 
pression, see Appendix [B]) Motivated by the observation that the formation rate of certain (though 


ees (For the derivation of this ex- 
not Earth-like) planets is proportional to the second power of host star metallicity , we consider 
the case m = 2, as well as m = 1) 

Another possible anthropic condition comes from the fact that if a halo is too dense, it may not 
host a habitable solar system because of the effects of close encounters [23]. Following Ref. [18], we 
assume this anthropic condition to take the form 


il 
My TVe ia (19) 


>There is observational data for metallicity of galaxies in our universe [22], which our crude model here does not 
reproduce quantitatively. However, when we straightforwardly extrapolate the empirical data to other universes, the 
same regions of the integrand in Eq. are suppressed /enhanced so that the effect on our calculation is qualitatively 
the same. It must be noted, though, that the strength of the effect may change; e.g. P(pa) for our model with m = 1 
is qualitatively similar to the distribution obtained using the observation-motivated method with m = 3. 


where n,, 0}, V;, and t,, are the density of stars, critical “kill” cross section, relative velocity of 


. ~ ; 1/6 
encounters, and some timescale relevant for the condition. Since ny x Pyir, Vy ~ Uvir «K M 1/ 3pi/6 


and o; and t,, are (expected to be) independent of M and p,;,, this is translated into 


n(Mt) = © (Bans — puelr(M)) (57) ). (20) 


where O() is the step function (= 1 for x > 0 and = 0 for x < 0), and we have normalized M by 
M ie: 
The value of Pmax is highly uncertain. One way to estimate it is to follow Ref. and take 


= Prir 2 Tvemw f NE? par \ 78 
nem (1pe)°(PA), ay ards thw te af EO ( MY (Boe) ay 
Re) Pvir,. Mw i ! Mp \Muw Pvir,. Mw ey 


where m, is the proton mass, ray © 1.5 x 108 km is the Sun-Earth distance, and pyirmw ~ 
2x 10-*° g/cm’, Tyirww ~ 5 x 10° K, and Myw ~ 1 x 10'°Mo are the virial density, virial 
temperature, and mass of the Milky Way galaxy, respectively. Using ta ~ tevo) = 5 Gyr, Eq. (19) 


leads to 
Muw 


2/7 
Bani S10 priest ( ) ~3 x 10-2 g/cm’. (22) 


This corresponds to the constraint from direct encounters, i.e. the orbit of a planet being disrupted 
by the passage of a nearby star. There can also be a constraint from indirect encounters: a pass- 
ing star perturbs an Oort cloud in the outer part of the solar system, triggering a lethal comet 
impact [18]. For a fixed M, this constraint can be about four orders of magnitude stronger than 
Eq. 

Pmax ~ 3 x 10~*° g/cm; (23) 


namely, our Milky Way galaxy may lie at the edge of allowed parameter space. 

In our analysis below, we consider either or both of the above conditions Eqs. (18) and (20). In 
the real world, there are (almost certainly) more conditions needed for intelligent life to develop. 
However, incorporating these conditions would likely improve the prediction/postdiction for pa. 
In this sense, our analysis may be viewed as a “conservative” assessment for the success of the 
framework, although it is still subject to uncertainties coming from the modeling of observers. 


5 Distribution of the Cosmological Constant 


Our modeling of observers has several parameters which need to be determined phenomenologically: 
Eq. (£6) contains Minin, Mmax; tevol, aNd thum, while Eq. (20) contains pmax. We take the “minimum” 
galaxy mass appearing in Eq. (16) to be 


Mmin = 2 x 10'4Mo, (24) 


P(pa) P(pa) 


=== 


We PA,obs 


Figure 1: The normalized probability distribution of the vacuum energy P(pq) as a function of pa/pa obs: 
The left panel shows P(p,) with the metallicity condition, Eq. (18), with m = 0 (i.e. no condition; dashed, 
blue), m = 1 (dot-dashed, red), and m = 2 (solid, black). The right panel shows P(p,) with the upper 
bound fmax, Eq. 20), with pmax = oo (i.e. no constraint; dashed, blue), 6 x 10776 g/cm? (dotted, purple), 
4.5 x 1076 g/cm? (dot-dashed, red), and 3 x 107° g/cm? (solid, black). 


below which the efficiency of star formation drops abruptly [17]. For teyo1, and tpumn, we take them 
approximately to be the age of the Earth and lifetime of the Sun, respectively: 


tevol = 5 Gyr, thumn = 10 Gyr. (25) 


In our analysis below, we use Eqs. (24) and (25); we do not impose the constraint from galaxy 
cooling, i.e. we set Mmax = 00. While the values of these parameters are highly uncertain, our 
results are not very sensitive to these values. The dependence of our results on them will be 
discussed at the end of this section. 

In Fig. [1] we present the normalized probability distribution for the vacuum energy P(p,) as a 
function of p/PA,obs, Under several assumptions about the function n: 


(i) “minimal” anthropic condition: n(M,t) = 1 
(ii) metallicity condition: Eq. with m = 1 and 2 


(iii) maximum virial density condition: Eq. With Pras ={3X10-"", 45x10"? 6x10") ¢/ em", 
which are {1,1.5,2} times the value in Eq. (23). 


(The result with pmax given by Eq. is virtually identical to the case with the minimal anthropic 
condition.) The left panel presents the effects of metallicity, showing (i) and (ii), while the right 
panel those of Pmax, with (i) and (iii). 

Interestingly, in all cases, our predictions prefer a positive cosmological constant over a negative 
one, as opposed to the situation in earlier measure proposals where strong preferences to negative 
values have been found [IJ]. In Table[I] we provide the probabilities of having p, > 0 (and < 0) in all 
six anthropic scenarios. The absence of an unwanted preference towards negative p, is satisfactory, 


No condition 
Metallicity, m = 1 
Metallicity, m = 2 


poise = 6% 10-* g/om® 
Praise = 245% 10-" g/ ame 
Paxec=3 *10-* gem 


Table 1: The probability of observing a positive and negative cosmological constant, P(p, > 0) and 
P(p, <0), for six assumptions on the anthropic condition. In all cases, a positive value is preferred over 
a negative one, consistent with observation. 


P(pr)PalPayoos P(Pr)PA/Pabs 
0.5 0.8 


0.6; 


087 1 10 100 Todo PAs 


1 00 1000 PA,obs 


Figure 2: Same as Fig. [1] but the horizontal axis now in logarithmic scale. To show the probability density 
per tenfold, the vertical axis is chosen to be pa P(pa)/Pa obs. The distributions are normalized in the region 
pa > 0. 


especially given that the measure of Ref. was not devised to cure this problem. It comes from 
the fact that the present measure does not have a large volume effect associated with the global 
geometry of anti-de Sitter space, which was responsible for a strong preference for negative p, in 
earlier, geometric cutoff measures [11]. In contrast with these measures, the quantum measure 
of Ref. [5] does not count the number of events; rather, it gives quantum mechanical weights for 
“situations,” i.e. quantum mechanical states as described from the viewpoint of a single observer 
(geodesic). The preference towards a positive value comes from the fact that for py > 0 some 
observers still form after vacuum energy domination, while for py < 0 it is not possible due to the 
big crunch. 

Figure [1] shows that P(p,) is always peaked near p, = 0, with the distribution becoming wider 
as the anthropic condition gets weaker. In Fig. 2] we plot the same distributions in logarithmic 
scale for pa/Par ops, limiting ourselves to py > 0. To show the probability density per tenfold, the 
vertical axis is chosen as p, P(p,)/Par ps. From these figures, we find that our anthropic assumptions 


10 


P(pa) P(pa) 
0.15 0.15; 


tevol 0.10+ tburn 
0.05 
20 30 at Pe ag 20 30 a Pe 


Figure 3: The normalized probability distribution P(p,) with the metallicity condition, Eq. (18), with 
m = 2. In the upper-left panel, Min is varied as 2 x 10'! Mg (solid, black), 6 x 10!‘ Me (dot-dashed, red), 
and 0.67 x 10''M. (dashed, blue); and in the upper-right, Mmax as oo (solid, black), 10'4Mo (dashed, 
blue), 10! Mo (dot-dashed, red), and 2 x 10!*Mz (dotted, purple). The lower left and right panels vary 
tevol and thurn as {(solid, black), (dot-dashed, red), (dashed, blue)} = {5,1,0} Gyr and {10,7,15} Gyr, 
respectively. 


lead to results that are consistent with the observed value within one or two orders of magnitude. 
In particular, metallicity alone is enough to bring the agreement to an order of magnitude level. 
This is because mergers, which lead to an increase in metallicity, are suppressed for larger values of 
pa due to earlier vacuum energy domination. This result is comfortable, especially given that the 
constraint from encounters is effective only if fmax is close to the Milky Way value, as in Eq. (28). 
Given our crude treatment of observers, we consider these results quite successful. 

Finally, we discuss the sensitivity of our results to variations of Mmin, Mmax, tevol, and thurn, 
which can be thought of as “systematic effects” of our analysis. In Fig. [3] we show the distributions 
of P(pa) with the m = 2 metallicity constraint, varying the values of Mmin, Mmax, tevol, aNd thumn, 
respectively. We find that, while the detailed shape of P(p,) does change, our main conclusions are 
robust: (i) There is no strong preference to a negative vacuum energy; in fact, a positive value is 
preferred. (ii) The predicted distribution of p, is consistent with the observed value at an order of 
magnitude level with the metallicity constraint. 
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Figure 4: The normalized probability distribution P(p,) with a metallicity condition: Eq. with m = 2. 
The light and dark shaded regions indicate those between 1 and 20, and outside 20, respectively. The 
observed value pa/PA,obs = 1 (denoted by a vertical line) is consistent with the distribution at the lo level. 


6 Conclusions 


In this paper, we have studied the probability distribution of the cosmological constant (or the 
vacuum energy pa) in the multiverse, using the quantum measure proposed in Ref. [5]. We have 
found that this measure does not lead to a strong preference for negative pa, as opposed to earlier 
measures proposed based on geometric cutoffs, because it does not experience a large volume effect 
associated with the global geometry of anti-de Sitter space. Moreover, we have found that a positive 
value of pa is preferred, consistent with observation. 

We have found that a simple, intuitive condition based on metallicity is enough to reproduce the 
observed value of pa, at an order of magnitude level. This is comfortable because effects from other 
possible anthropic conditions, such as the ones from encounters, are much more sensitive to the 
details of the conditions. In Fig. [4] we present the normalized distribution P(p,) with the m = 2 
metallicity constraint, where the 1 and 20 regions are indicated. We find that the observed value is 
consistent with the calculated distribution at the 1a level. 

It would be interesting to refine our analysis including more detailed anthropic effects, such as 
those of star formation. Another possible extension of the analysis is to vary other cosmological 
parameters, such as the primordial density contrast Q and spatial curvature Q, (at a specified time), 
in addition to px. We plan to study these issues in the future. 
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A Press-Schechter Formalism and Fitting Functions 


The Press-Schechter function F’ is 


where 6,(t) and o(M,t) are given by 


1.629 + 0.057 e~2:3Gweat for pa > 0 
d-(t ~ 2.5 11 
1.686 +0.165(4-) +0.149(—4-) for pa <0, 


torunc 
with tounch defined in Eq. (£7), and 
o(M,t) ~ Qs(M) Git). 
Here, Q is the primordial density contrast, 
s(M) ~ (9.1 ie ae {50.5 logy9(834 + pt) — 92°" : : 

where ps = Meg? with £ = Pmeaties/ iy 2.7 SV, and 

0.20655 jtanh?/* (SH qt){1 —tanh'?"(3 Hat) on 

PA 
G(t) ~ + 1.437{1 — cosh“/9(3Hat)} | for px > 0 
4/3(41/342/3 | 2 oie 
0.549 £4/3GN342/ [1 +0.37 ( ) ja - ( ) for pa <0, 


terunch terunch 
where Hy = \/87Gy|pa|/3. 


The function H in the extended Press-Schechter formalism is given by 


I, _ MME AD pte. sg / 
H(t; M,t) =— ot ant M8) dM’, 
— 1 Se(ts)  60(to) 
_ c\4l _ Ye 2 
B( My, ti, Mo, te) = ee TCA ERTAD) faa ae)), 


with s(M) and G(t) defined in Eqs. and (80). 

The virial density as a function of time can be fit, following Refs. [25], as the density evolution 
of a closed universe, according to Birkhoff’s theorem. The virial density is then given in terms of 
the density at turn-around rescaled by the ratio of the volumes, pyir = (Ryir/Rturn)?Ptum- Here, 
Ryir/Riun — 2 at early times (t < 1/H,) as well as for |p,| — 0 at any fixed t. For positive 
Pa, Rvir/ Rium = 2/(V3 — 1) & 2.73 at late times 26], while for negative pa, Ryir/Rium — 27/° for 
t > teunch- Our fit is given by 


ah. 

9 sinh? (3 Hat) 1.41 ‘oli 1.41 

(187 Pmaiter (0) Tee” ) + (40.8 Pa) , for PA > 0 (33) 
Pvir a 2/3 

2 sin =) 123.6 

(18x Pmarter() era 123.647 (e* oe a 


4 rerumck —1) 
where Pmatter is the matter energy density. This fit is accurate to better than © 5% and 2% for 
pa > 0 and < 0, respectively. 
Finally, the time at which most of galaxies of mass M forms, i.e. the solution to Eq. (15), is well 
approximated by the following fitting function: 
1.880 + c(a) M + 3(a) M3 +¢5(a) M®> for —105 a. 


3/2 ~ a ,obs 
T(M)/Gyr ~ = co(a) + c4(a@) M + c4(a) M* + c(a) M* for 0 = A < 10 (34) 
cy(a) +c{(a) M + c¥(a) M3 + ch(a) M° for 10 < az S 100, 


where a = (Pay Pnzorn) QO one Q)*, M= logio SANS and 


c(x) = —0.311 + 1.276 82?" + 1.412 log. {1 + |x|%"}, 
c3(t) = 0.470 — 0.656 e%78* — 0.317 log, (0.2 + |z\), 

c5(t) = —0.0142 + 0.0381 €%7* + 0.00822 log, (0.05 + |z]), 
ch(x) = 1.880 — 0.00205 2, 


i} 
o(x) 
‘(z) = 0.408 + 0.569 e~ 1 + 0.295 log, (1 + 2), 

c(z) = 0.277—0.251e~* — 0.125 log,,(1+ 2), (35) 
L(x) = —0.000889 + 0.0151 e~* — 0.00220 log,9(1 + 2), 


cl(x) = 1.880 — 0.00205 x, 

c'(x) = 0.767 — 0.00293 x — 230.274, 

ci(z) = —0.530 + 0.000336 x + 0.847 2-1, 

c!(z) = 0.106 — 0.0000118 x — 0.125 x~°! — 0.0131 log,)(—5 + 2). 


°For pa < 0, the approximation leading to Eq. (83), i-c. pyir 2 a few (Pmatter+pa), breaks down for t/terunch & 0.8, 
where we should rather use pyir = 0 (since there is no stable structure forming). However, since pyir in Eq. is 
small there anyway, using it up to t/terunch = 1 does not lead to a significant error. 
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This fit is accurate to better than ~ 5% for M & 10!!M,z (but it becomes worse for smaller M, e.¢., 
the accuracy is + 12% at M ~ 6 x 10!°Mo). For 100 < pa/paops < 150, we use the last expression 
of Eq. (84), good to the level of ~ 10%; and for pa/Paobs > 150, we use 


pvir(T(M)) /(10776 g/cm*) ~ (3.66 + 0.032 a) — (1.36 + 0.0013 a)M, (36) 


which is accurate to the level of © 10% up to pa/Paobs © 4500. 


B- Anthropic Condition from Metallicity 


In this appendix, we derive the function n arising from the metallicity constraint, Eq. (I8). Suppose 
that in a merging tree of a galaxy 7 at time t, 7 is found to have progenitor galaxies 7 = 1, 2,--- 
with varying masses Vj; at time t’ < t. Note that this also includes accretion, i.e. matter that was 
not part of galaxies of appreciable size, since F(M = 0,t) = 1 in the Press-Schechter formalism, 
where accretion is treated as mergers of extremely tiny galaxies with a large galaxy. 


Now, let us assume that the relative mass fraction in j7 that came from 7 and 7’ is given by 
dF(M,,t!) ) dF (My!) 

dM; / dM) 
respectively. This is true within the Press-Schechter formalism as long as M;v < M,, since then 


, Le. the ratio of total amount of baryons at time ?¢’ in galaxies of type 7 and 7’, 


the overdensities within spherical top-hat regions containing masses M;,;, and M; are independent 
of each other at early times. Once Mj, ~ M;, the assumption is not justified, but in these regimes, 
there can only be a small amount of merging occurring from 7,7’ to 7, implying little contribution 
to metallicity. The assumption, therefore, provides a good approximation. 

Let x; (i = 1,2,---) be the fraction of baryons in the universe that formed stars in halos of 
mass M; at time t’. In our simple model, the star formation rate is proportional to the rate of halo 
formation for masses M > Myin and otherwise zero: dx;/dt' « O(M; — Mnin) d?F(M,,t’)/dM,dt'. 
The increase in total metal content summed over galaxies of mass M; is taken to be proportional 
to the star formation rate therein, dx;/dt, so the increase in (linear) metallicity dZ; is 


dF(M,,t')\~* da; at)\ @ int! 
( 1) dx: ys (sa) TPMEt) out, — Moin) dt’, (37) 


dZ,(t') «x dt' « 
e) ( dM; ae dM; dM,dt' 
i.e. the total increase in metal content divided by the total mass. The increase in metallicity of 
galaxy j due to stars at time t’, then, has to be weighted by the relative matter fraction of galaxies 
i, as described above: 

dF(Mi,t') dF (Mi,t!) 
Di a EZi("’) ye dM, dt! O(Mi — Mmin) ait (38) 


> dF(M;,t") ms > dF(M;,t) D 
i dM; i dM, 


dZ; (t’) = 


where we must normalize to the total mass of galaxy j at each time t’. 
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In the continuum limit, the sum over 7 becomes an integral over masses. Therefore, the metal- 


licity of galaxy 7 of mass M at time t is 


F pM yy py d2F(M' 2’) i / / 
2(M,t) o | ag to tM Sara O(M" ~ Manin) _ fe el F (Monin!) — FOE} (a9) 
0 fovam' ae!) 0 1— F(M,t’) 


where ¢ = min{t, top}. The constraint that one cannot accumulate negative metallicity determines 
the timescale a as a solution to 


d 


AF (Minin, t') — F(M,t')} = 0, (40) 


thats 
at which time merging of galaxies of mass Minin < M’ < M into those more massive than M begins to 
dominate over formation of new galaxies in this mass region. Since merging into larger structures is 
also occurring at earlier times, one expects that we slightly underestimate the metallicity. However, 
in practice, the formation of new galaxies in this mass range and mergers into structures beyond 
are well separated in time, so a simple cutoff at Laie is sufficient. 

Now, since F(M,t') < F(M, tstop) = erfe(1/V2) ~ 0.317, the denominator of Eq. is always 
between 0.68 and 1; in fact, it is very close to 1 in most of the parameter regions. Therefore, we 
can safely ignore the denominator of Eq. and obtain 


Z(M,t) x (F(Mnin; #) — F(M,#)). (41) 


In general, the probability of forming planets is expected to be proportional to some power m of 
the metallicity [21]. This gives 


n(M,t) = Z(M,t — tev)”, (42) 


which is Eq. in the text. 
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